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ABSTRACT 



The standard paradigm of hierarchical structure formation in a ACDM universe predicts the presence of dark matter subhaloes, hosted 
by Milky Way-sized galaxies. Anticipated subhalo masses range from 10'" down to a cut-off mass between 10"^ and 10"" Mq. If 
dark matter is composed of heavy self-annihilating or decaying particles, these subhaloes could be visible in the y-ray band as faint 
and temporally constant sources without astrophysical counterparts. Based upon realistic subhalo models and current observational 
constraints on annihilating dark matter scenarios, we predict that one massive Galactic subhalo between 10^ and 10^ M© may already 
be present in the 11-month catalogue of Fenni-hAT. Indeed, at least twelve objects in the first Fermi catalogue qualify as candi- 
dates. The most promising object, IFGL J0030. 7+0724, is investigated in detail using a dedicated Swift X-ray follow-up observation 
and a refined positional analysis of the 24-month Ferini-hAT data. With the new observations, seven point-like X-ray sources have 
been discovered, of which SWIFT J0031 19.8-1-072454, which coincides with a faint radio source (12mJy at 1.4 GHz), serves as a 
counterpart candidate of IFGL J0030.7+0724. The broad-band spectral energy distribution is consistent with a high-energy-peaked 
blazar. However, flux and extent of IFGL J0030. 7+0724 may also be compatible with a dark matter subhalo. Detection of temporal 
variability or improved astrometry of IFGL J0030.7+0724 are necessary to rule out or confirm an astrophysical origin. We discuss 
strategies to identify y-ray sources that are associated with self-annihilating dark matter subhaloes. 
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1. Introduction 

Several astrophysical observations indicate that in the early as 
well as in the present Universe a non-baryonic form of dark mat- 
ter (DM) prevails over the baryonic matter content. Structure 
formation favours a cold dark matter (CDM) scenario (for re- 
cent reviews see, e.g., iBertone et al]l2005H D^ Amico et a ni2009t 
iBertondbOlOh . However, the nature of DM remains unknown. 
A class of promising candidates for CDM are stable, weakly 
interacting, massive particles (WIMPs) with masses between 
10 and 10^ GeV, predicted by theories that extend the standard 
model of particle physics. The most prominent extensions en- 
compass those based on supersymmetry and universal extradi- 
mensions, which were invented to solve inconsistencies of the 
standard model at high energy scales [C>(TeV)], and which de- 
liver adequate DM candidates in this way. These particles can 
self-annihilate or decay, producing detectable signatures in the 
final states such as energetic photons (y rays), antimatter, and 
leptons. 

Unravelling the nature of DM remains a challenging 
problem for astronomy and particle physics, and a variety of 
attempts to detect signals have been made, using both direct and 
indirect detection techniques. For instance, multi-wavelength 
observations of astrophysical targets have constrained the 
self-annihilation rate of DM, which is related to the thermally 
averaged annihilation cross section. In particu lar, regions with 
high D M densities such a s the G al actic Centre (lAharonian et al] 
I2006allbl iMorselli et all I20T0I: lAbramowski et alj 1201 IblT 



Galactic Ridge (lAharonian et al. 2006dl). dw ar f spheroidal 



galaxie s (dSphs; 



2008b; Aharonian et al 
Aharonian et al. .2010; 



._ „ j Albert et al.l 

'200 81 120091; lAliu et alJ 120091: 



Colafrancesco et al 
T 



l2007t 



Abdo etal.1 l2010cl; lAcciari et all 



20 101; lAleksiaet all 1201 iL IH.E.S.S. Coll aboratio n 201 lb. "as 



well as globular (IColafrancesco et alJ 2006c iWag ne 
Abramowski et al.| 201 lal) and galaxy clusters jAleksic et d] 



20 1 Ot lAckermann et all 1201 Oh 



serve as excellent targets. 
Furthermore, DM annihilation in the entire Galactic halo as 
well as its subhalo population produces a diffuse y-ray flux, 
which contributes to the overall diffuse signal of the Galaxy. 
Comparatively stringent upper limits on the annihilation cross 
section have been obtained from the combination of both dedi- 
cated observations and boun d s obtained from the diffuse y-ray 
flux (lAbazaiian et alJ 120101; ICirelli et al.l l2010t iMeade et all 



l2010l:|Papucci & Strumiall2010l:IZahariias et alj|2010l) 
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Based on the theory of hierarchical structure formation, DM 
haloes of Milky Way-sized galaxies are anticipated to host nu- 
merous DM subhaloes with m asses between a cut-off scale 
10 " -10-^^ and lO'^M^ (e.g.. lBringmannll2009l) . where Mq 
denotes the solar mass. This expectation is a consequence of 
the early collapse of overdensities in the expanding Universe 
(Diemand et al. 200l), leading to the formation of initially low- 
mass haloes, which subsequently serve as building-blocks for 
larger haloes b y merging at later times. Besides analytical cal- 
culations (e.g.. iBerezinskv et al.ll200l l2006l l2008l) . recent nu- 
merical high-resolution A^-body simulations of stru cture forma- 
tion in a ACDM cosmology (Komatsu et aT]l20in) . such as the 
Aquarius P roject ( Springel e t al. 200 8a||biror the Via Lactea II 
simulation dDiemand et al. 2008bi ; iZemp et alJl2009h . allow to 
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study substructures in detail. For a Milky Way-type galaxy, these 
simulations predict a large number of subhaloes (up to lO'^) 
with masses distributed following a power law, AN I AM oc M^", 
where a e [1.9; 2.0]. The DM density profiles of large-mass sub- 
haloes are found to be similar to the ho st's, which results in high 
central densities jSpringel et al.ll20"08ah . The spatial distribution 
of subhaloes is "anti -biased", i.e., the dominant fraction is placed 
far away from the host-halo's centre. 

In self-annihilating DM scenarios, subhaloes are expected 
to appear as weak point-like or moderately extended y-ray 
sources, and a small fraction of them could be detectable 



scopes (e.2..|Pieri et a 


.1120051 120081 1201 It Snringel et al. 


12008b; 


Kuhlenetal.1 120081 


Andol 120091 iBucklev & Hoopei 


201C; 


Brun et al.l 201 ll IZechlin et al.l 201 1 ). Examples for currently 



operating telescopes are the Fermi-LAT (20MeV - 300 GeV, 
lAtwood et al.l |2009|) and imaging air Chere nkov telescopes 



(lACT s; E > l OOGeV) such as H.E.S.S. dAharonian et al.l 
l2006d) . M AGIC (lAlbert et all l2008at iTridon et al.l 1201 Oh . and 



VERITAS ("Wee kes et al.ri2002h . In the near future, a signif- 
icant improvement in the overall sensitivity and lower en- 
ergy threshold will be achieved by upcom ing experiments 
such as H.E.S.S. -II (Iyincentll2005l) and CTA (lHermannll20Tol: 



iThe CTA Consortiumll2010t lDoroH201 iK Such instruments are 



possibly sufficiently sensitive to detect nearby large-mass sub- 
haloes of (5(10^) Mo within distances of 0(1) kpc. 

This paper consists of two separate parts, which can in prin- 
ciple be read independently. The first part (Sects. |2]to|4]i inves- 
tigates the detectability of subhaloes with Fenni-LAI, where the 
basic framework for predicting the y-ray properties of subhaloes 
is laid out in Sects. |2| and [3] In Sect. |4] properties of detectable 
subhaloes are investigated by means of a fiducial source. In the 
second part. Sect. |5] discusses the search for DM subhaloes in 
the first Fenni-LAI point-source catalogue (IFGL) and subse- 
quent multi-wavelength studies of the most promising candi- 
date, IFGL J0030. 7-1-0724. A discussion of the physical origin 
of IFGL J0030.7-H0724 and prospects for lACTs are presented 
in Sect.|6| 

Throughout this paper, Hubble's constant is Hq - 
73kms 'Mpc yielding the present value of the Universe's 
critical density p„it = 3//,^/(87rGN) - 1.48 x 10"MoMpc"\ 
where Gn denote s Newton's gravitational constant (e.g., 
ISpergel et al.ll2007l) . 



2. Gamma rays from DM subhaloes 



On the other hand, a second and more realistic model is consid- 
ered to account for subhalo evolution, henceforth referred to as 
subhalo model (SHM). 

2.1. Density profile 

The subhalo's DM density profile p(r) is assumed to follow 



p(r) 



Ps 



{rlr,)y{l+rlr,f 



1 for r < Tcut, 
for r > rcut, 



(1) 



In the following, the y-ray flux from DM subhaloes will be de- 
rived, based upon current theoretical models of the correspond- 
ing radial density distribution. 

With respect to undisturbed, isolated galactic haloes, hence- 
forth field haloes, the general formation history of (embedded) 
subhaloes differs significantly. Analytical models and numerical 
A^-body simulations of structure formation found their physical ln(c™'^) = 
properties to depend on particular evoluti onary conditions, i.e., ~^ 
formation time, evolution, and orbit (see iDiemand et ani2007l 
[2008a, and references therein). Tidal interaction with the grav- 
itational potential of the host halo leads to tidal stripping and 
heating, and can therefore truncate the outer region of subhaloes. 
In the following, two different approaches will be discussed. On 
the one hand, subhaloes are modelled assuming negligible tidal 
effects and are therefore considered to be in a genuine virialised 
state. Because this approximation is (at least) valid for field 
haloes, this model will be tagged as field-halo model (FHM). 



where r denotes the distance to the subhalo's centre. In gen- 
eral, the profile cuts at an outer radius rcut, which is the virial 
or tidal radius (/?vii or Rx), respectively. Given y = 1.0 for 
the remainder, the profile follows the universal spherically sym- 
metric Navarro-Frenk- White (NFW) profile, well-fitting h aloes 
resolved in numerical simulations'! (iNavarro et al ] ll997h . The 
profile is defined by two parameters: a characteristic inner ra- 
dius Ts, where the effective logarithmic slope of the profile is 
-2, and an inner density ps - A-pir^). In case of FHM haloes, 
which are not subject to tides, both parameters are related to 
each other by the virial halo mass Myir. This quantity is defined 
as the mass inside the sphere of radius 7?vir, which encloses a 
mean density of Ac times t he critical density o f the Universe 
at the considered redshift z ( iNavarro et al. I IT997I IB ullock et al.l 
l200ll) . Mvii := 47r/3 AcPciit^Jij.. The virial overdensity at z = is 
Ac ~ 100, as suggest ed by models of t he dissipationless spher- 
ical top-hat collapse (lEke et al.lll996l: iBrvan & Normanlll998l) 
and assuming present concordance cosmology. In general, the 
subhalo mass M is given by a volume integration of Eq. [Tj re- 
vealing M = Anp^rlf{c), where /(c) = ln(l + c) - c/(l + c) 
and c denotes the concentration parameter of the subhalo. For 
non-disturbed haloes, the concentration is then given by the 
virial concentration = /?vir/'"s- Generally, the concentra- 
tion depends on the subhalo mass and redshift, c = c(M, z), 
where lighter haloes have higher concentrations ( INavarro et al.l 
,1996, 1997; Bullock et al. 2001) . Since observation al estimates 
are lacking (see Sect. 2.2.1 in iLavalle et al.l I2O O8'. and refer- 
ences therein), c{M) is ad opted from A^-body simulations. For 
the FHM, the toy model of iBuUock et al.1 (1200 llF I is used, where 
the halo's (average) virial concentration at redshift z is connected 
with the density of the Universe at the halo's (mass-dependent) 
collapse redshift zdM), Cvk- = K{1 + Zc)/{1 + z). The contrac- 
tion parameter K is constant and independent of cosmology. 
To determine the mass dependence of Zc and, therefore, of Cyir 
at z - 0, the low-m ass extrapolation of the Bullock model by 
iLavalle et al.l (l2008l) is adapted. Because Cyh implicitly depends 
on Ac (see also Sect. I2.2l i. a con version of Cyjr to Ac = 100 
was appliec0 with the relation of iHu & Kravtsovl ( llool . The 
concentration-to-mass relation is well-fitted by the polynomial 
form 



FHM 



In 



M 



(2) 



cf™ = {4.265, -0.0384, -3.91 x 10"^ -2.2 x 10"^ -5.5 x 10"^). 
Note that this model almost equals the relation derived by 



Note that details on the very inner slope of halo profi les remain to 
be clarified, by simulations as well as observationally (e.g.. lWalker et al.l 

2011:Sa lucci et aL.2007.') . 

^ This model extends a proposal by INavarro et al.l ( 119971) . 
' ILavalle et al.l chose = 81.6. 
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Fig. 1. Concentration-to-mass relation for subhaloes at z = 0. 
The dashed blue line depicts the concentration predicted by 
the Bullock model (FHM), while its corresponding scatter is 
given by the blue-shaded area. The mean concentration cor- 
rected for subhalo evolution (SHM) is shown by the solid red 
line for a galactocentric distance of 74k pc, the average dis- 
tance of the dSph galaxies included from lAbdo et al. I (120103) 
(black triangles). For comparison, the concentration derived for 
Aquarius subhaloes is indicated by the dot-dashed dark green 
line within its validity range, the low/high-mass extrapolation 
by the double-dotted green line. The black filled circle marks 
the virial concentration of the MW. 



iPieri et aT for a cosmology as used in the Aquarius sim- 

ulation. Regarding the concentration of SHM haloes, the low- 
mass extrapo lation of the Bu llock model provides a conservative 
estimate (cf.. iPieri et al.ii[2008i) . 

However, because subhalo formation differs from that of 
field haloes and includes tidal truncation at Rt, the vi rial concen- 
tration is not well defined for subhaloefl ( iDiemand et al. 200^. 
Therefore, the SHM incorporates an empirical correction of Cvir- 
Indicated by numerical simulations, the concentration of sub- 
haloes increases with decreasing distance to the host's centre 



:) 



El 

^ VIT 



(3) 



see IDiemand et al.1 (I2007ll2008ah and lKuhlen eTaP (l2008h . The 
galactocentric distance is Dgc = (R^ + — 2RoD cos / cos b)^^^, 
where D denotes the subhalo's distance to the Sun, (/, b) its posi- 
tion in galactic coordinates, an d j^n = (8.28 ±0.2 9) kpc the Sun's 
distance to the Galactic Centre dCatena & Ulliol2 010). The virial 
radius of the Milky Way (MW) is R^J^ = c^* rf * ^ (288 + 
61) kpc, where c!; ^^ ^ 19.70+2.92 and rf * = (14.65+2.24) kpc 
dCatena & UmdMloh . The power-law slope ao - 0.237 is 
adopted as fitting subhaloes resolved in the Aquarius simulation 
dPieri et al I I20TI . 

Intrinsic to the stochastic process of halo formation, the con- 
centration of individual haloes scatters around the median c pro- 



^ In general, the physical subhalo radius R, is smaller than the for- 
mally defined virial radius i?vir, implying the physical subhalo mass M, 
to be smaller than My,^. For massive subhaloes, the M,(Mvir) relation is 
approximately linear, where M,/Mvi, ~ 0.2 (see AppendixlAt. Given the 
empirical model correction discussed below, the formal virial quantities 
will be used in the remainder of the paper. 



vided by the quantities c™^(M) and c™'^(M, Dgc), respectively. 
The corresponding probability distribution follows a lognormal. 



P(c,c) 



log 10 e 



V2^< 



■ exp 



1 ^logio 



(4) 



where (Ti^^^^, = 0.14 dBullock et al.ll200lHWechsler et al.ll2002h . 

The concentration-to-mass relations are shown in Fig.[I] For 
the FHM, the scatter is also depicted (68% c.l. of logjo c), see Eq. 
|4] In addition to the concentration of the MW, Fig. [T] contains 
a selection of eight dSphs that are associated with sufficiently 
precise stellar data, which allow a con clusive modelling of the 
DM distribution (see lAbdo et al.ll2010d) . Each dSph is modelled 
with a NEW profile with parameters chosen to fit measurements 
of stellar line-of-sigh t velocities and their distributions (see also 
iMartinez et al.ll2009l) . The dSph's virial concentration is given 
by its characteristic density, ps = AcPcritcJ;j./[3/(Cvir)], where 
tidal effects on the inner system are assumed to be negligible. 
The SHM is depicted for D„c = 74 kpc, the average galacto- 
centric distance of the dSph subset. Additionally, the models are 
confronted with direct predictions of the Aquarius simulation, 
derived from scaling relations fitting subhaloes observed in the 
simulation. Details are provided in Appendix lAl 

Within its scatter, the concentration model of FHM haloes 
consistently describes the DM profile of dSph galaxies and 
the MW itself. However, the median values c™^ underpredict 

' vir ^ 

dSphs, whereas the subhalo model SHM provides convincing 
agreement (as expected by A'-body simulations). The concentra- 
tion derived directly from the Aquarius simulation confirms the 
SHM within the validity range, see Fig. [1] Note that the mean 
distance of subhaloes resolved in Aquarius is 64 kpc. 

2.2. DM annihilation in sublialoes 

For self-annihilating particles, the total rate of photons (or par- 
ticles) emitted by a DM subhalo with energy E in the interval 
[Ei;E2] is 



{o-v}effNy 

2mi 



dVp^(r) oc 



El jm^ 



r dNy 



(5) 



£1 /m^ 



where (o-v)e(f is the thermally averaged annihilation cross sec- 
tion times the relative velocity, the WIMP mass, and dNy/dx, 
X = Ejm^, denotes the differential spectrum of photons per an- 
nihilation. Assuming «: D, the produced photon flux is given 
hy 4> - -Lli^TiD^). The solution of the integral holds for y = 1 .0 
and c » 1 . In Eq. |5] a small, flat core replacing the u nphysical 
singularity at the halo centre dBerezinskv et al lll992h is safely 
neglected (given the NEW profile used here). For y - 1.0, Eq. |5] 
simplifies via = [3M/(47rAcPcritcJ:J] 



1/3. 



7(Cv.r)2 



18m2 



(6) 



For ay = 1.2 profile (Eq. [T]i the photon rate increases 
by a factor of ~ 1.5 for subhaloes above IO-'Mq. DM an- 
nihilation in subhaloes may be additionally boosted by sub- 
subst r ucture populations (s ee IStrigari et al.l 120071 iKuhlen et al.l 
I2OO8I: IMartinez et alj|2009t) . Conveniently, die value of <crv)eff 
is normalised to the value (crv)o = 3 x lO^^^cm^s"', which 
leads to the correct relic density. An increase of the annihila- 
tion rate, a so-called boost factor (crv)eff/(crv)o, could in princi- 
ple be related to the underlying particle physics framework (e.g. 
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bb, 500 GeV 
W*W", 500 GeV 
tV, 500 GeV 
tV, 150 GeV 



100 



10 




E [GeV] 

Fig. 2. Differential y-ray energy spectra E"^ ^ originating from 
final-state fragmentation of WIMP annihilation. Four different 
WIMP models are depicted: mass = 500 GeV, final states: 
bb (solid black line), (dashed red line), t^t" (dotted blue 

line); mass - 150 GeV, final state t^t" (dot-dashed green 
line). The parametrisations are valid down to E/m^ ^ 0.01. 
The grey-shaded area indicates the considered energy range 
Ee [10; 100] GeV. 



iFornengo et al. 20041) and effects such as Sommerfeld enh ance- 
ment (e.g.. lArkani-Hamed et alJl2009l;lKuhlen et alJl2009h . 

We considered DM to be composed of self-annihilating 
heavy WIMPs of mass = 500 GeV and investigated three 
distinct annihilation models: two of then]_with total annihila- 
tion in heavy quarks or gauge bosons {bb and W^W") and a 
model with total annihilation in the leptons r^r . Additionally, 
we considered WIMPs of = 150 GeV for annihilation in 
T^T" final states. Particles of this type are, for instance, pro- 
vided by supersymmetric theories, e.g., manifested in the neu- 
tralino. The WIMP masses chosen are compatible with WIMPs 
which might explain the recen tly observed c osmic -ray electron 
and positron excess, see, e.g., iMeade et alj ( |20IO|) . In general, 
heavy WIMPs are also supported by collider searches such as 
the non-detection of supersymmetric particles in the 7 TeV run 
of the La rge Hadron Col l ider (for an integrated lum inosity of 
35 pb-i) jAad et al.ll20TTl; ICMS Collaboration! I20T1I and refer- 
ences therein). 

Given these final annihilation states, hadronisation and the 
subsequent decay of ;r"-mesons lead to a continuous y-ray spec- 
trum. The resulting photon sp ectra dNy/dx were mo delled using 
parametrisations provided bv lFornengo et alJ (l2004l) . see Fig. |2] 
Note that photons produced by final state radiation (FSR) and 
virtual internal bremsst rahlung (VIB) (.Bergstrom et al...2005a<ib| ; 
iBringmann et alJl2008l) are neglected, because a significant con- 
tribution of FSR is only expected for W^W" at high energies 
(E > 0.6 m^) and possible contributions of VIB are highly 
model-dependent. 



3. Candidate gamma-ray sources 

Given a high WIMP mass, a DM subhalo will show up as steady 
(very) high-energy y-ray source. The differential photon spec- 
trum follows a hard power law (index F < 1.5) that cuts off 
exponentially at energies beyond 10 GeV, see Fig.|2l 

Candidate sources are selected according to their possibil- 
ity to originate from DM subhaloes based on their observational 



quantities flux and angular extent. Note that the detailed spectral 
shape of faint sources is observationally rather unconstrained. 
Via Eq.|6] the effective self-annihilation cross section (crv)eff re- 
quired to obtain a given flux (p for the intrinsic source extent 0^ is 
determined hy X.- A-nD^ip, where 9^ constrains the distance D to 
the subhalo. For feasible candidate sources the required (crv)eff 
should comply with current observational constraints. 

Conveniently, the characteristic profile radius (see Eq. [1]) 
traces the intrinsic extent of a DM subhalo, because for an NFW 
profile 87.5% of the total luminosity is produced within r^, (see 
Table [1] for values of r^). Therefore, the distance to the subhalo 
is D a: r^jO^, where 9^ denotes the angle corresponding to r^. 
Owing to the centrally peaked profile, about 68% of the total lu- 
minosity is emitted within the angle 068 - 0.46 9s. The following 
relations are given with respect to 0s and can easily be adjusted 
for 068, which is more convenient for a comparison with obser- 
vational data. With - [3M/(47rAcPcritcJij.)]'^'', the distance to a 
subhalo with given 0s is related to its mass and concentration. In 
the FHM, 



Ofhm(M;0s)-3.8 



M 



lO^Mp 



1/3 



/ „FHM \ 



37.9 



kpc. 

deg/ 



(7) 



Note that the concentration depends on the subhalo mass via 
Eq.|2]as well as, in the SHM, on the object's position (/, b). 

For a given WIMP model (crv)etf is then fully determined by 
the subhalo mass (Eq. |6l) and the observed quantities flux and 
extent: 



<crv)eif(M;</.,0s) = 96 7r-- 



Ny \4AcPc: 



5/3 



cp M-"^f(c,irf 
0? c'. 



(8) 



Additional contributions to the DM signal from annihilation in 
the smooth halo as well as the entire subhalo population were 
neglectecQ The required (crv)eff is highly sensitive to the (obser- 
vationally unconstrained) concentration, because Eq.|8] roughly 
depends on c^r'^.. 

4. Interpretation of Fermi sources as DM subhaloes 



Based on the study of a fiducial candidate source in Sect. 14. II the 
properties of Fermi-hPS detectable subhaloes are investigated in 
Sect.lO 



4.1. A fiducial candidate 
4.1.1. Observational properties 

In com bination with improv ing (integrated) sensitivity at high 
energy dAtwood et al J 12009^, and cf.. Fig. |7]), the expected en- 
ergy spectrum of DM subhaloes (Sect. |2] and O favours a de- 
tection at the high-energy band of Fenni-hPil. In Appendix 
B, we investigate the detection sensitivity for faint, moderately 
extended (0s < 1°, corresponding to 068 ^ 0.5°), and high- 
latitude {\b\ > 20°) sources between 10 and 100 GeV in detail. 
We find that a spectrally hard high-energy source with a flux 
1^(1 0-1 00 GeV) = 1.4xl0"'°cm ^s ' and moderate extenfl 0s = 
1° can be detected as a point-source with a reconstructed flux 



^ For the fiducial candidate in Sect. |4] this additional contribution is 
less than 1%. 

* In comparison with the point spread function of Fermi- 
LAT, 6s = 1° (as implying 068 ~ 0.5°) corresponds to about 
3o"psp, where o"psf ~ 0.15° for energies beyond 10 GeV (see 
http://www-glast.slac.stanford.edu/software/TS/glast_lat_performance.htm). 
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M [M , J M [M : ] 



Fig. 3. Effective self-annihilation cross section (crv)efr required for a moderately extended, faint Fermi-LAT source to originate from 
a DM subhalo of mass M. Assumed source parameters are 0(lO-lOOGeV) = 1.6 x lO"'" cm^^s"', an intrinsic extent 6^ = 1° 
(068 = 0.5°), and the Galactic position {l,b) - (114°, -55°). The dashed blue and solid red lines indicate the average prediction 
considering the FHM and SHM, respectively. The intrinsic SHM scatter is shown by the red-shaded area and the dotted red lines. 
The left panel depicts WlMPs of = 500 GeV totally annihilating in bb, while th e right panel consi ders niy = 150 GeV wit h total 
annihilation in t^t . Current contraints on {crv)^s from lPapucci & Strumial(l2010l) (grey-shaded) and lAbazaiian et all (1201 Oh (Ught 
grey-shaded) are plotted in combination with the expectation from thermal freeze-out (blue-shaded). 



0p(lO-lOOGeV) = 0.9 x lO^'^cm^^ ^-i^ with a sky-survey ex- 
posure of one year. With respect to the true flux (p emitted by the 
entire source, the reconstructed flux (^p fitted by the point-source 
analysis in general decreases with increasing 9^. To account for 
this effect, the scaling relation (p(0^) - h{9^)(p^ is used in Eq.|8] 
where h(9^) = 1 for 9, « 2 crpsp and h{9^) * 0.72 (ejdeg) + 0.89 
for extended sources up to ~ 1° (see AppendixiBlfor details). 

The high-energy flux of the fiducial source above 10 GeV 
has been chosen to be at the level of the detection sensitiv- 
ity, 0p''( 10 - 100 GeV) = lO^'^cm^^s"', assuming an extent 

of 9^^ - 1°. Given the dependence of the SHM concentration 
on the galactocentric distance (see Eq.O, the fiducial source is 
placed on a particular line-of-sight chosen to match the location 
of IFGL J0030.7+0724 (investigated in Sect.|5l). In general, this 
line-of-sight serves as an appropriate (conservative) benchmark, 
because it points to (anticentric) positions where the majority of 
subhaloes is located. 



4.1.2. Subhalo interpretation 

Adopting the properties of the fiducial source. Fig. [3] depicts the 
effective enhancement factors {crv)^ffl {crv)^ required to generate 
the emission (fi^^ by DM annihilation (obtained via Eq.[8]l. In the 
left panel, WlMPs of - 500 GeV are considered to totally 
annihilate in bb, while the right panel assumes - 150 GeV 
and annihilation in r^r . For a given WIMP model, the resulting 
enhancement factors of the FHM and SHM are widely different. 
With respect to the FHM, much less enhancement is required in 
the SHM, which is manifested in generically higher concentra- 
tions of SHM subhaloes. Within the scatter of the concentration 
intrinsic to the stochastic nature of halo formation (Eq.lUi, which 
is shown for the SHM, the necessary enhancement spans about 
one order of magnitude. Only moderate enhancement is required 
for massive subhaloes between 10^ and 10^ M©, where the low- 
est {crv)efF is needed for - 150 GeV and t^t" final states 
(amongst the WIMP models considered here). The lowest pos- 
sible enhancement factors within the concentration scatter of a 



10^ Mq subhalo are listed in Table[T]for the different subhalo and 
the WIMP models of Sect. |2] 

The distance to the fiducial candidate anticipated in the FHM 
and the SHM is shown in Fig.|4] The intrinsic concentration scat- 
ter implies a corresponding distance scatter for a given halo mass 
and angular extent. Note that a similar scatter is present for the 
FHM, but is not shown in the figure. Compared with the FHM, 
tidal effects lead to higher concentrated subhaloes. This in turn 
favours a closer distance at the same mass and angular extent for 
SHM subhaloes than for FHM. 

4.1.3. Consistency witli observational constraints 

The resulting values of (crv)eff can now be checked for con- 
sistency wi t h the d iffuse extragalactic y-ray background (EGB; 
lAbdoetal.ir2010bh . Fermi-LAT measurements of the overall 
diffuse y-ray flux allow the derivation of the isotropic high- 
energy EGB, which is shown to be compatible with a feature- 
less power-law spectrum (T - 2.41 ± 0.05) and integrated dif- 
fuse flux (pECB(> 100 Me V) = (1.03 ± 0.17) x lO^^ cm-^s-'sr-' 



Table 1. Enhancement factors (crv')etf/(crv)o required to explain 
the fiducial y-ray source with a DM subhalo of 10^ Mq. 



Model 


= 500 GeV 

bb W+W- T^T- 


150 GeV 

T+T- 


[kpc] 
















FHM 


> 321 > 479 


> 1386 


> 138 


0.067 


y 37.86 j 


SHM" 


> 8 > 12 


> 35 


> 3 


0.029 


j ^SHM \ - 1 
1 86l6 ) 



Notes. The factors correspond to the WIMP models discussed in the 
text. The respective photon yields are A'*'''"''*"'^'^'(10 - 100 GeV) = 
6.95 (4.66) [1.46] for m^. = 500 GeV and ^^"(10- 100 GeV) = 1.61 
for nijf = 150 GeV. We list the minimum values within the c-scatter. 
In addition, the subhalo's (average) characteristic radius rs(M) = 
rl[MI(\Qf' Mo)]'/V--(cvi,) is depicted. The normahsation of cj™ im- 
plies the Galactic position (114°, -55°) and intrinsic extent 6^ = 1°, 
corresponding to £) » 1.7 kpc (ZJgc ~ 8.8 kpc). 
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Fig. 4. Distance to the fiducial subhalo in the FHM (dashed blue 
line) and SHM (solid red line), respectively, as function of the 
subhalo mass M. The scatter of the SHM distance is indicated 
by the red-shaded area. An extent of 0s = 1° (^68 ~ 0.5°) and a 
source position (I, b) - (1 14°, -55°) are assumed. 

Table 2. Upper limits on (cr\')eff /(crv)o from the EGB. 



WIMP model 


Upper limit on {crv)^^ 


/<o-v>o 


Channel 


[GeV] 


FHM 


SHM 


FHM 


SHM 


bb 


500 


1650 


530 


1875 


605 


w+w- 


500 


2096 


673 


2381 


769 




500 


3490 


1121 


3964 


1279 




150 


378 


121 


429 


139 




10- 


"'Mo 


10 


-'Mo 



Notes. We assumed a subhalo mass fraction of /sh = 15 % for a cut-ofF 
mass M,^i„ = 10"' Mq. Upper limits are list ed for the cut-ofF m asses bor- 
dering a 500 GeV neutralino scenario (see iBringmannI [20091) and with 
respect to lA = 180° and E = 40 GeV. See text and Appendix [Cl for 
details. 



dAbdo et al.ll2010bh . The diffuse y-ray flux anticipated from DM 
annihilation in the Galactic halo as well as the entire subhalo 
population is shown in Fig. |5] (see Appendix ICl for details), in 
comparison with the EGB. Both FHM and SHM subhaloes are 
depicted for - 150 GeV and annihilation in t^t . In this 
model, the nearly isotropic diffuse flux from the subhalo popula- 
tion contributes about 1% to the EGB (assuming SHM subhaloes 
and no sub-substructure) and is fainter than the contribution of 
the smooth halo (> 3%). Note that the flux from the extragalactic 
halo population is lower than the contr ibution of Galactic sub- 
haloes, see, e.g.. lAbazaiian et 30(120101) . 

The contribution from the smooth halo peaks at the Galactic 
Centre, where a high astrophysical foreground is also present, 
and can therefore not be isotropic. Given that the EGB has 
been derived assuming isotropy, the most robust upper limits on 
(crv)eff are determined by the subhalo contribution and are listed 
in Table |2l depending on the WIMP model and cut-off mass. 
The bounds were obtained requiring that the specific intensity of 
the subhalo population (7^(180°, £)) does not exceed the EGB, 
where {/v(iA, E)) depends on the angle kjj between the Galactic 
Centre direction and line-of-sight and the y-ray energy E (see 
AppendixO. 

However, more stringent constraints have been provided by a 
more detailed modelling of the EGB, including all DM compo- 
nents. To evali iate a possible DM orig in of the fiducial source, 
the results of lAbazaiian et al.l (1201 Ol) and iPapucci & Strumial 




80 100 
V [deg] 

Fig. 5. Average diffuse y-ray flux from DM annihilation in the 
Galaxy as function of the angle between the observational and 
Galactic Centre direction i^. The fluxes from the host halo 
(double-dot-dashed magenta line), FHM (dashed blue line), and 
SHM subhaloes (solid red line) are depicted. Each curve was 
derived at the peak energy of vly(E) (40 GeV) assuming total 
annihilation in t^t" with - 150 GeV and {crv)ss - {crv)o. 
The minimum subhalo mass used was M,nin - 10"^ M© (see 
Appendix0. The EGB at 40 GeV is shown by the dotted black 
line. 



(l2OT0lfl are included in Fig. [3] As stated in the introduction, com- 
petitive and similar constraints have been also provided by the 
non-detection of various objects with high (central) DM densi- 
ties. Note, for instance, that the constraints u sed here are consis - 
tent with recent bounds from dSph galaxies dAbdo et al ] l2010ch . 

4.2. Fermi-LAT detectable subhaloes 
4.2.1. Expected number 

Integrating over the mass and spatial distribution (Eq. IC.2! ) re- 
veals the total number of detectable subhaloes with masses Mi e 
[M; M + AM], solar distances Di € [D;D + AD], concentrations 
Cl 6 [c;c + Ac], galactic latitudes bi € [b,b + Ab], and galactic 
longitudes / 6 [0; 27r], N = NiMi, Di, ci, bi), 



N = JdM J iDD-J dc j ib cosb J ilP{c, 



c) 



d»sh(£>gc, M) 
dM 



(9) 



Mi Di 



Cl bi 



In 



^aNfdMAlfde,^^^-^ f dc c"^ fdfo cos fdZ f (c, c) — , 
J J sin^e, J J J ' ' ' dM' 

Mi 'c„,m(e„M) bi 

using D = rs/tan0s (see Sect. |3]l and = 3/(47rAcPcrit)- 
Parameters defining the subhalo distribution were taken to match 
the Aquarius simulation (Appendix 0. In total, this resulted 
in about 6.4 x 10'^ Galactic subhaloes residing in the Galaxy. 
For every single 6^ and M, the integral counts detectable sub- 
haloes only, i.e., their concentration is sufficiently high to en- 
sure their (crv)eff to be smaller than the observational constraints 



^ In compa rison with IPapucci & Strumi^ the work by 
lAbazaiian et al.l includes a fore- a nd background su btraction. Note that 
the MW halo parameters used by lAbazaiian et alj are similar to those 
adopted in this work (see Sect. l2.lt . 
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Fig. 6. Expected number of one-year detectable subhaloes per 
mass decade at \b\ > 20°. The subhaloes' mass and spatial dis- 
tributions were adopted from Aquarius (Appendix ICll while the 
SHM concentration was used. Mmin - 10"^ M©. The number 
corresponding to the = 500 GeV (t^t", - ISOGeV) 
WIMP model is shown in the left (right) panel. The mass as- 
signed to each bar denotes the geometric mean of the interval. 
Red and red-patterned bars show t he number cons i dering ob- 
servational constraints on {crv)eff bv lAbazaiian et al] ( 1201 Ol) and 
iPapucci & Strumial (1201 Oh . respectively. 



(cf.. Fig. [3]). Therefore, the lower bound of the concentration in- 
tegral Cmin(ds,M) is determined via Eg. El choosing the instru- 
ment's sensitivity and constraints from lAbazaiian et alj (1201 Ol) 
and Papucci & Strumia (201(|), respectively. To account for the 
fact that highly extended objects will be hardly detectable (see 
Appendix |B]), we conservatively chose 6*™"" = 1°. For the 
SHM, Fig. |6] shows the number of detectable subhaloes per 
mass decade expected in one year of data taking while consid- 
ering subhaloes at galactic latitudes \b\ > 20° only (cf.. Sect. 
|5]|. Given the dependence of Cmin on the WIMP model, the re- 
sults for the bb, - 500 GeV model are compared to the 
niy - ISOGeV scenario. 



Considering the constraints found bv lAbazaiian et al ] (120101) 
for the bb (t^t") model, on average 0.4 (0.2) subhaloes within 
10^ and 10^ Mq are anticipated for detection with Ferm/-LAT 
in one year. Given the Poisson distribution of A^, this means 
that up to one massive subhalo is expected in the one-year data 
set (at 95% confidence ). Comparable res u lts have been c laime d 



by other authors, e.g., 
LPierietal. (2008, 2011 



Baltz et al.1 (|2008|). iKuhlen e t al.1 (l2008h . 



, and lAnderson et a l. (2010). 
In general, note that numerical simulations like Aquarius 
and Via Lactea II neglect the influences of bary onic matter dis- 
tribute d in galactic disks. A recent study by iD'Onghia et al.l 
(I2010h indicates that a baryonic disk may reduce the number of 
(massive) subhaloes in the inner galaxy by a factor of 2 to 3. 

4.2.2. Properties 

Given a particular y-ray source (such as the fiducial of Sect. 14. 11 1, 
massive subhaloes between about lO*" and 10** M© require a min- 
imally enhanced annihilation cross section (crv)eff, see Fig. |3] 
Consistently, the probability for these objects to appear in cur- 
rent data sets peaks for high subhalo masses (Fig.|6]l. Therefore, 
subhalos with masses of 10^ up to 10** M© in corresponding dis- 
tances from O.S to 10 kpc (Fig.|4]i are favoured for detection with 



Fermi-LAJ' as faint and moderately extended sources. Although 
for the FHM a DM origin of the fiducial source is excluded, 
within the scatter of the more realistic SHM cross sections re- 
quired for sources with (^(10- 100 GeV) * 10 '°cm"^s ' and 
angular extents up to ~ 1° are well consistent with observa- 
tional constraints. The presence of sub-substructure will even 
lower {crv) ,ff by a m ass-dependent fac tor of ~ 2 to 3 for massive 
subhaloes dKuhlen et al. 2008; Marti nez et al I2OO9I) . In case of a 
cuspier profile (7=1 .2), the required cross section is lowered by 
an additional factor of 1 .S. At least for WIMPs of m;^, = ISO GeV 
annihilating to t*t", this leads to a required enhancement of the 
order of unity within the scatter Summarising, in optimistic but 
realistic scenarios a y-ray emitter at the (one year) detection level 
of Fermi-LKT with a measured extent 9^% of about O.S° might be 
consistent with a subhalo driven by self-annihilating DM. 

In the next years, Ferm;-LAT will provide deeper observa- 
tions with increased observation time Tobs. This will improve 
the sensitivity by roughly VT^, leading to a factor of about 2 
lower values of the minimum detectable flux for the five-year 
catalogue. This in turn will allow us to detect fainter subhaloes 
with a correspondingly reduced minimum {crv)^^. The average 
number of detectable subhaloes within five years is about 1.3 
(0.8) for the Z7^, = SOOGeV(T+T", = ISOGeV) scenario. 

Via Eq. [8j a comparison with observational constraints on 
(cv)eff allows to estimate the maximum flux subhalo candidates 
are expected to have. The catalogued flux is 0p oc 0^ / h{0^), which 
is quadratic for 0^ <K 0.3° and linear in the limit of large 0s- 
The increase with 6^ originates from decreasing subhalo dis- 
tance. Because massive subhaloes require minimum (crv)eff, a 
subhalo of 10^ M© is assumed below. Within the concentration 
scatter, the maximum flux between 10 and 100 GeV expected for 
a source with 6^ - 1° is 



^'"(10- 100 GeV) ^ 1.2(2.8)x 10"" 
^ <o-v>o 



11 (cn')eff _2 -1 
cm s 



(10) 



in the bb, niy - SOOGeV and t^t 



= ISOGeV sce- 



nario, respectiv ely. Given the observational constraints of 
lAbazaiian et al.L the high-energy flux of catalogued candidates 
should not exceed 0™" < 4.0 (2.8) x 10 



5. Searches for DM subhaloes in the 1FGL 

In the previous section we demonstrated that DM subhaloes 
could appear in y-ray catalogues of sufficient sensitivity as faint, 
non-variable, and moderately extended objects without astro- 
physical counterpa rts. The 1 1-montlB point-source catalogue of 
Ferm/-LAT (IFGL. lAbdo et al.ll20I0al) fists 14S1 sources signif- 
icantly detected above 100 MeV (test statistic TS > 2S, corre- 
sponding to a significance S - 4.1cr), together with the flux in 
five discrete energy bins (up to 100 GeV), position, significance 
of variability, and spectral curvature. Source spectra have been 
fitted with power laws. Among the sources, 630 object^ are not 
confidently associated with known sources at other wavelengths. 

Although sophisticated methods have been applied to find 
multi-wavelength associations for unidentified sources, all al- 
gorithms suffer from lacking sensitivity or incomplete sky- 
coverage of current surveys. Therefore, the sample of unasso- 
ciated high-latitude Fermi-LAT sources is expected to be com- 
posed of several source classes, among them faint AGN (Active 



August 2008 to July 2009 
' Adding sole associations with other y-ray catalogues, 671 sources 
are "unassociated". 
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Table 3. Cut efficiencies on the sample of unassociated sources 
and AGN. 



Cut 


Unassociated 


AGN 




671 


693 


\b\ > 20° 


249(100%) 


539 (100%) 


non-variable 


241 (97%) 


372 (69%) 


detected between 


12 (5%) 


58(11%) 


10-lOOGeV 



Notes. The cuts are cumulative, i.e., each number includes all cuts listed 
by previous rows. See text for details. 



Galactic Nucle i), galaxy clusters, and new exotic sou rces like 
DM subhaloes (IStephen et al.l2O10HMirabal et al.l20To[ and ref- 
erences therein). Concerning the IFGL catalogue, i mproved as- 
sociation methods recently presented by Mirabal et al.l (l2010t) 
revealed that < 20% of all unassociated IFGL sources with 
\b\ > 15° may contain new y-ray emitters. 

To single out possible subha lo candidates within the sam- 
ple of unassociated sources (cf., iB ucklev & Hoopeill2010h , we 
searched the sample for non-variabl43 sources detected between 
10 and 100 GeV. Requiring a detection at high energy provides 
subhalo candidates driven by heavy WIMPs and avoids confu- 
sion with high-energy pulsar^ Furthermore, the candidate's lo- 
cation was constrained to galactic latitudes \b\ > 20° to avoid a 
general confusion with Galactic sources. Additionally, the lower 
Galactic foreground i mproves the detectio n sensitivity of Fermi- 
LAT at high latitudes dAtwood et al.ll2009h . 

Applying all cuts, twelve unassociated sources remain. 
The twelve sources are listed in Table |4] together with addi- 
tional information from the catalogue. With the exception of 
IFGL J0614. 1-3328, the sample consists of sourc es at the faint 
end of the entire IFGL sample. Given the result of lMirabal et al.l 
(1201 Ol) . the sample should statistically contain two to three sub- 
haloes at most, consistent with the estimate discussed previously 
(Fig. |6]l. The expectation of the sample consisting mostly of 
AGN is met by applying the same cuts to all AGN detected 
by Fermi-LAiT. A comparison with the sample of unassociated 
sources reveals similar cut efficiencies (5% vs. 11%, see Table 
O, indicating that the two populations share common proper- 
ties. Note that for the AGN the variability cut has subdominant 
influence as well, see Table [3] Except for three, all AGN that 
passed the cuts have been classified as BL Lac. 

Even though the twelve candidate objects are listed in the 
IFGL catalogue as unassociated, we extended the counterpart 
search to a wider choice of astronomical catalogues. Table|4]lists 
the classifications of counterpart candidates in the 68% confi- 
dence regions around the most likely IFGL positions, retrieved 
from the NASA/IPAC Extragalactic Database (NED). In par- 
ticular, possible associatons are provided by radio and X-ray 
sources, since most of the selected y-ray sources are expected 
to be AGN. Given that no detailed association study was con- 
ducted, some of the tabulated sources might be by-chance asso- 
ciations. 

Governed by lacking association, faintness, and spectral 
shape, this study focusses on the most promising candidate, 
IFGL J0030.7 4-0724. Within the errors, its high-energy flux and 

The cut is passed by sources with a steadiness probability Ps > 1%. 
" The spectral pr operties of y-ray pulsars can mimic the spectra of 
DM subhaloes, see lBaltz et al.l ( |2007|) . Ho wever, spectral cut -off ener- 
gies of y-ray pulsars are well below 5 GeV jAbdo et aU2009iJ) . exclud- 
ing a detection above 10 GeV. 



sources and candidates in the first Fermi-hAT catalogue 

10-8 [ — ^ 

IFGL J0030.7+0724 




0.1 1 10 100 



E [GeV] 

Fig. 7. Energy spectrum of IFGL J0030. 7+0724 as given in 
the first Fermi-hPS catalogue between 100 Me V and 100 GeV 
(filled circles). The solid red line depicts the (catalogued) power 
law fitting the data. Only the highest energy bin has been signif- 
icantly detected, while for the low-energy bins upper limits are 
shown (95% c.l.). 



spectral index are well-compatible with a self-annihilating DM 
scenario. The source has only been detected between 10 and 
100 GeV, see Fig. 13 

5.1. Multi-wavelength properties of 1FGL J0030.7+0724 

5.1.1. Catalogued data 

No counterpart candidate was found within the positional 
uncertainty of the y-ray source at a 68% confidence level 
(Table I?}. In the 95% confidence region, the faint radio object 
NVSS J003030-I-07 2132 is located (/i.4GHz = (3.5 ± 0.4) mJy; 
Condon et al. ''1998'). However, no conclusive infrared (2MASS, 
Skrutskie et al. 2006) or optical (USNO Bl.O, Monet et al. 2003|) 
association of the NVSS source is known so faQ Note that 
no dSph galax y is located in the source region (NED). ROSAT 
(0.1-2.4keV: IVoges et al.ll999h observations of the region with 
an exposure of about 170 s revealed no X-ray source down to an 
energy-flux level of ~ 10^'^ergcm^^s"' (Borm2010). 

We emphasise that the orphaned faint radio source is likely 
located in the uncertainty region of IFGL J0030.7+0724 by 
chance, because about 0.7 NVSS sources are expected by statis- 
tics. 

5.1.2. Ferm/-LAT data 

By analysing the 24-month public archival data between 10 and 
100 GeV, updated results on IFGL J0030.7+0724 will be pro- 
vided. For the same energy range, a reanalysis of the 1 1 -month 
data is presented for comparison. Particular focus will be drawn 
on positional properties, the high-energy flux, and the photon 
distribution, which allows us to investigate possible counter- 
parts, temporal variability, and the angular extent. 



'2 Within the 2cr positional uncertainty of NVSS J003030-I-072132, a 
very faint optical SDSS (Sloan Digital Sky Survey) s ource is located 
(26.0"') - SDSS J00303 1.22+072 132.2 (SDSS DR7. lAbazaiian et al.1 
200§). However, this object was observed with the edge of the plated 
SDSS camera. Therefore, this detection is probably spurious. 
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Table 4. Unassociated, non-variable IFGL sources at high galactic latitudes. 



Name 


0"68/cr95 


S 


/p(0. 1-100 GeV) 


r 


0p(lO-lOOGeV) 


Possible associations" 


Remarks 


IFGL J 


[arcmin] 


[0-] 


[10"" erg cm"- s"'] 




[10-"'cm-2 s-'] 


68% c.l. 




0022.2-1850 


6.0/9.6 


9.4 


1.3(4) 


1.6(1) 


1.6(7) 


RadioSs (4,21,22), Gs (20 - 18) 




0030.7-H0724 


3.0/5.1 


5.8 


1.0(4) 


1.7(4) 


1.5(7) 


— 




0051.4-6242 


2.4/4.2 


12.0 


1.8(5) 


1.7(1) 


1.7(8) 


Gs (20), XrayS (3.8) 


c 


0143.9-5845 


3.0/4.7 


9.0 


1.4(4) 


2.0(2) 


2.0(9) 


Radios (28^), Gs (20- 13) 




0335.5-4501 


2.4/4.0 


8.6 


1.5(4) 


2.1(2) 


1.6(8) 


Gs (19,18) 




0614.1-3328 


1.2/1.7 


54.4 


11.2(6) 


1.93(3) 


3(1) 


GrayS 


b 


0848.6+0504 


5.4/8.6 


5.4 


1.0(5) 


1.2(3) 


1.6(8) 


RadioSs (2,3,5), Gs & *s, XrayS (4.4) 


c 


1323.1+2942 


1.8/2.7 


11.9 


1.5(4) 


2.0(1) 


2.1(8) 


RadioSs (2.8,263,724), Gs & *s 




1754.3+3212 


2.4/4.1 


15.6 


2.6(4) 


2.09(9) 


1.4(7) 


Radios (38+) 




2134.5-2130 


3.0/5.1 


6.7 


1.1(3) 


1.9(2) 


1.4(7) 


Radios (22), Gs (20) 




2146.6-1345 


3.0/4.4 


9.8 


1.5(5) 


1.8(2) 


1.8(8) 


Radios (23), Gs (20), XrayS (1.9) 


c 


2329.2+3755 


1.2/1.9 


10.4 


1.7(5) 


1.6(2) 


2.4(9) 


G(14) 


c 



Notes. The columns list the positional uncertainty (Tf,8(95) [68% (95%) c.l., semimajor axis], detection significance 5 in Gaussian sigma, integrated 
energy flux /p(0. 1 - 100 GeV), spectral index T, and the photon flux (/ip(10-100 GeV). Here, parentheses indicate the corresponding error on the last 
decimal(s). Furthermore, the type classifications of sources found in astronomical catalogues within the 68% uncertainty region of the Fermi-hKI 
position are listed. Classifications referred to are RadioS (radio source), G (galaxy), * (star), XrayS (X-ray source), and GrayS (y-ray source). 
For radio, optical, and X-ray sources corresponding fluxes are given in mJy (at 1.4GHz ['^*: 843 MHz, 4.85 GHz]), apparent magnitudes, 
and 10"'- ergcm"^ s"', respectively. The unabsorbed X-ray flux was derived from the catalogued count-rate, assuming a power law with index 
2.0 (with WebPIMMS, http://heasarc.gsfc.nasa.gov/Tools /w3pimms.html ). T he hydrogen col umn density was obt ained from the LA B survey, see 
Sect l5X3l Sources referre d to are listed in the FIRST JSecker et al.)ll995!), JVAS/CLASS jjackson et al1[2007l). NVSS , SUMSS jMauch et al.l 
l2003h . 2MASS, APMUKS jMaddox et al.ll990l) . SDSS ([Abazaiian et al. 2009), ROSAT, or EGRET jHartman et all 199^ catalogue, respectively. 

The spectrum is probably curved. The y-ray so urce has been associa ted by a cross-correlation of unidentified Fenni-hAT sources with the 
ROSAT All Sky Survey Bright Source Catalogue (see lStephen et alj | |2010l) for details). 



Table 5. Positional and spectral properties of 
IFGL J0030.7+0724 as given in the catalogue (11 months) 
and by the (re-)analysis of the first 1 1 and 24-month data sets. 



Data 


E 


R.A. 


Dec. 


<^p(10-100GeV) 


5(0p) 


set 


[GeV] 


(J2000) 


(J2000) 


[lO-'^cm-'s-'] 


[T] 


IFGL 


0.1-100 


00 30 42.6 


+07 24 09 


1.5 ±0.7 


6.6 


11 


10-100 


00 30 37.6 


+07 24 15 


1.4 ±0.7 


6.5 


24 


10-100 


00 30 47.6 


+07 24 20 


0.9 ± 0.4 


6.6 



Notes. The second column lists the analysed energy range. The 1 1 (24)- 
month analysis focusses on the high-energy flux (/ip(10- 100 GeV) only. 
In all cases, the significance S of the high-energy bin is well above 6a. 



The data analysis was performed with the latest public ver- 
sion of the Fermi Science Tools (v9rl8p6l3 along with recom- 
mended options and the set of instrument-response functions 
P6.V3J)IFFUSE (Rando 2009). Throughout the analysis, the 
optimiser MINUIT was used. For reliable results, photons of 
event class 3 (Diffuse) and 4 (DataClean) within a radius of 10°, 
centred on the nominal position of IFGL J0030. 7+0724, were 
selected. Given that most of the photons are at highest energies 
(Fig.|7]i, only events between 10 and 100 GeV were selected to 
minimise the background and ensure a narrow PSF The data 
were processed using gtselect, gtmktime, gtltcube, gtexpmap, 
and unbinned gtlike. To compute the most likely position and 
its corresponding uncertainty on basis of the 10- 100 GeV 
photon sample, we used gtfindsrc. For the purpose of detailed 
counterpart searches, the two-dimensional likelihood function 
L(R.A.,Dec.) was computed, which provides the 95% uncer- 
tainty contour by 2 A(logL) = 6.18 (2 degrees of freedom). 

The source model for the data analysis contains all 
IFGL sources within the region of interest (ROI, ra- 

Fermi Science Support Center, |http://fermi.gsfc.nasa.gov/ss(71 



Table 6. High-energy photons detected from 
IFGL J0030.7+0724 within 0.5°. 



E 


R.A. 


Dec. 


1? 


Af 




event 


[GeV] 


[deg] 


[deg] 


[deg] 


[30 d] 


CT 


class 


83.8 


7.6330 


7.3975 


56.26 


2.46 


B 


3 


11.8 


7.7293 


7.3771 


36.39 


5.19 


F 


4 


39.8 


7.7841 


7.4962 


47.38 


7.96 


B 


4 


10.2 


7.6426 


7.4483 


34.21 


10.46 


F 


4 


15.0 


7.6361 


7.1872 


38.24 


11.12 


B 


4 


43.8 


7.8392 


7.4151 


20.81 


18.93 


F 


4 



Notes. The table lists their energy E, celestial position (J2000), incli- 
nation )?, detection time Af, and conversion type (CT). By Af the time 
between detection and mission start is given. The conversion type is 
front (F) or back (B). For each event, we list the classification assigned 
by LAT data reconstr uction (Pass 6), wh ere 3 tags the Diffuse and 4 the 
DataClean class (see lAbdo et al.ll2010bh . 



dius 10°). Their parameters were taken as catalogued 
and we used the latest Galactic {gllJem_v02.fit) and ex- 
tragalactic (isotropic Jem_v02.txf) diffuse background mod- 
els. All parameters but those of IFGL J0030.7+0724 were 
kept fixed. Furthermore, the catalogued power-law index of 
IFGL J0030.7+0724 was used while fitting the flux between 
10 and 100 GeV. Although the exposure of the 24-month data 
has almost doubled with respect to the catalogue, the use of 
the catalogued properties for sources within the ROI will not 
affect the analysis between 10 and 100 GeV: The three nearby 
sources, i.e., IFGL J0022.5+0607, IFGL J0030.4+0451, and 
IFGL J0023. 5+0930, are not only more than 2° away from 
IFGL J0030.7+0724, but they are also not significantly detected 
between 10 and 100 GeV. Furthermore, visual inspection does 
not reveal any other relevant source within this nearby region. 
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The analysis of the 1 1 -month data reproduces the catalogued 
values well (Table |5]i. After 11 (24) months, five (six) photons 
between 10 and 100 GeV have been detected within a radius 
of 0.5° around the nominal position, listed in Table |6] Except 
one, all photons are classified as class 4 events and are therefore 
very likely signal events. The Galactic foreground and the ex- 
tragalactic background at the source position are negligible with 
respect to the signal, with an expected total number of back- 
ground photons A^bg = 0.6 (1.2) within the considered region of 
0.5°. For comparison, the predicted number of signal events is 
Nsig - 4.9 (5.8) after 1 1 (24) months. According to the 1 1 -month 
data set, the (10-lOOGeV) best-fit position shifts by about 2.5'. 
The small positional error of the sixth photon also accounts for 
the increase of the source's positional uncertainty, see Fig. [8] 

The average flux over the entire data set has decreased by a 
factor of roughly 1.5 with respect to the first 11 months (Table 
EJ. To judge on the variability of IFGL J0030.7+0724, its tem- 
poral photon distribution (Table |6]l was tested for compatibility 
with a co nstant flux, using an unbinned Kolmogorov-Smirnov 
(KS) test dPress et al.ll2007h . The KS test is already vaHd for 
low photon counts, unlike the binned chi-square method used by 
the catalogue. The KS test confirms the null-hypothesis of non- 
variability with a probability of about 0.7 (0.5) for the 1 1 (24) 
months data set. The varying exposure on the region was taken 
into account by examining the photon distribution of the bright 
pulsar nearby (IFGL J0030.4+0451). 

The analysis of the (intrinsic) spatial extent of the source 
is based upon a likelihood-ratio test, using all photons listed 
in Table |6] The corresponding statistical measure is given by 
L(0^) = -2Y,t^ln[pd^t{Xi -x;0s) + b], where pdet{x;ds) is the 
probability distribution function for a photon detected at x, x 
denotes the best-fit position (Table |5]l, and b incorporates the 
flat background. For a spatially extended y-ray emitter /jjet - 
Ppsf * Pint, the (two dimensional) convolution of the Fermi- 
LAT PSF (P6_v3, diffuse class) with the intensity profile of the 
emitter In the subhalo case, the intensity profile follows the 
line-of-sight integral of the squared NFW profile (Eq. [TJ. The 
quantity AL — L - Lmin follows a chi-square distribution with 
one degree of freedom, with additional terms of the or der of 
l/jyi/z^ which are important for a small number of counts (IWilksl 
ll938HCashlll979l) . The likelihood is minimised (Ln,i„) for the in- 
trinsic extension parameter fitting the photon distribution best. 
Examining the 11 -month data, the test shows the source to be 
consistent with a point source, implying that the intrinsic extent 
is smaller than the (average) PSF (about 0.15°). The 24-month 
data favour a moderate extent Gs - O.H^jJ p deg, which is, how- 
ever, not significant. Upper limits on the extension parameter are 
0s 0.54 (0.72) deg at 95% confidence level, derived from the 
1 1 (24) months data set. Since the low statistics affect the chi- 
square d istribution, note that the confide nce level is not precisely 
defined (ICashlll979l:FMattox et al.ll 19961) . Furthermore, we point 
out that the PSF of Ferm;-LAT (P6_v3) may be underestimatecQ 
and changes will have an impact on the fitted extent. 
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Fig. 8. Celestial region of IFGL J0030.7-I-0724, 25' x 13'. The 
catalogued position is indicated by the red +, while the dashed 
red ellipse borders its uncertainty (95% c.l.). The black x marks 
the 24-month position, the solid black contour its uncertainty at 
95% confidence. Photons detected by Swiji-XKl (10. Iks) are 
mapped by the back-image, which is smoothed with a Gaussian 
(7"). The positions of NVSS radio sources are given by the dark- 
green boxes, the two NVSS sources discussed in the text are 
named. In this region, seven X-ray sources have been discov- 
ered, indicated by the blue arrows (see Table |2ll. Note that the 
boxes' size is chosen arbitrarily. 



in photon-counting (PC) mode. The XRT data were calibrated 
and selected with standard screening criteria {xrtpipeline), using 
the HEAsoft 6.10 package for data reduction with the current 
version of calibration files available (rele ase 2010-09-30). For 
the analysis, events with grades 0-12 (iBurrows et al ] l200l 
were us ed. The spectr al analysis was carried out with Xspec 
(12.6.0, lArnaudI Il996h , using the PC grade 0-12 response 
matrix swxpcOtol2s6J20070901v01 l.nnf with the ancilliary 
response function generated by xrtmkarf for PSF correction 
and the position of the source considered. The on-source region 
was selected to contain about 90% of the PSF (a: 47"). For 
background subtraction, an off-source region with radius of 
about 4' was used. To ensure a spectral fit of suflicient quality, 
the spectra were rebinned to a minimum of 10 events per bin 
(with grpphd). Owing to the low statistics accumulated, the 
C-statistic was used for spectral fitting. 

In the field-of-view (FoV) of XRT, seven new X-ray sources 
were discovered with a probability of being background fluctu- 
ations smaller than lO"*". We show them in Fig. [8] Their posi- 
tional properties, measured flux, and the flux corrected for pho- 
toelectric absorption between 0.2 and 2 keV are listed in Table 
Q The spectra of the two brightest sources are well-fitted by an 
absorption corrected power-law model, fixing the hydrogen col- 
umn density A^h to the nominal Galactic value. The power-law 
index for the faint sources was fixed to 2.0. The Galactic hy- 
drogen column dens ity was obtained from the LAB HI survey 
(iKalberla et al.ll2005l) for the corresponding celestial positions. 



5.1.3. Sw/ft-XRT data 

The field was successfully proposed for observation with 
die X-ray telescope (XRT, () .2 -10ke V) aboard the Swift 
satelhte dGehrels et alj |2004 l2005h . The observations 
(Obs. ID 00041265001) were carried out on 10 November, 
2010, between 00:23:46 and 19:52:56 UT with a total effective 
exposure of 10. 1 ks. Observations with the XRT were performed 



http://fermi.gsfc.nasa.gov/ssc/data/analysis/LAT_caveats.html 



Associations. The newly discovered X-ray sources were stud- 
ied for possible associations in other accessible wavelengths. 
Multi-wavelength surveys covering the region are the NVSS 
in the radio, the 2MASS in the infrared, and the USNO Bl.O 
and SDSS DR7 catalogues for the optical band. For every Swift 
source we found at least one SDSS source to be positionally 
coincident (Table |8]l, with apparent magnitudes between 21™ 
and 17™. Owing to insufficient sensitivity, the very faint SDSS 
sources have not been detected by USNO. 
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Table 7. X-ray sources detected with the Swift-XRT. 



ID 



Name 
SWIFT J 



[arcsec] 



S/N 



/'■""(O.l-lkeV) 
[10"'"* ergcm"^ s"' 



[10^" cm--] 



90 

[lO^^keV-' cm-2 s"'] 



[10- 



"'X0.2-2keV) 
erg cm"^ s"'] 



A 003000.3+072301" 6 3.5 

B 003017.8+072142 5 5.4 

C 003022.1+072623" 6 3.1 

D 003030.0+072013" 5 5.1 

E 003049.8+072316" 6 3.0 

F 003054.9+072328" 6 2.8 

G 003119.8+072454 5 6.5 




2.0 
1.4 ±0.3 

2.0 

2.0 

2.0 

2.0 
1.6 ±0.3 



5.2 ± 1.5 

6.7+2-3 

2 2+'-^ 

4 4+1:^ 

3 0+'-^ 

-U 
20 7+*^ 

-4 7 



Notes. The FoV is centred on (R.A., Dec.) = (7.6315, 7.421 1) deg with a radius of 13'. We give an internal ID, the position 
(SWIFT IHHMMSS.s±DDMMSS) and its coiTesponding error at 90% confidence level ergo (determined with xrtcentroid), and the signal-to- 
noise ratio SIN (Ximage) of the observed flux /"''^ If constraining, a power-law model corrected for photoelectric absorption was fitted to the 
spectrum, dtp/dE = (po (E/keWy'^. The hydrogen column density Nh was fixed during the fit. The unabsorbed flux f"""^^ was derived from the 
power-law fit. *"' Due to a low S/N a two-parameter power-law fit is not constraining. The fluxes were derived assuming the index F = 2.0. 



Table 8. Likely counterparts of the X-ray sources listed in Table|7] 



ID 


NVSS*"' 




2MASS 




USNOBl.O'"* 


SDSS DR7 






Name 


5 [miy] 


Name 


K 


Name 


R 


Name 


g 


Type 


A 


003000+072255 


47(2) 






0973-0005315 


20.0"' 


1003000.24+072254.7 


20.3"' 


6 


B 










0973-0005428 


20.4"" 


J003017.75+072140.6 


19.6"' 


6 


C 














J003022.22+072621.3 


21.4'" 


6 


D 






00302977+0720101 


15.3"" 


0973-0005481 
0973-0005484 


18.9"' 
19.3"' 


J003029.77+072010.3 


18.5"" 


3 


E 














J003049.61+072313.5 


21.0"" 


6 


F 






00305500+0723233 


15.7'" 


0973-0005560 


18.2"' 


J003054.80+072323.1 
J003055.00+072323.2 


20.7'" 
18.5'" 


6 
6 


G 


003119+072456 


11.6(6) 






0974-0005617 


18.6"' 


J0031 19.71+072453.5 


17.4'" 


6 



Notes. Scans ranging from radio (NVSS) to infrared (2MASS) and optical (USNO and SDSS) wavelength bands are given. The table lists the 
object's name and the catalogued flux or apparent magnitude. Here, parentheses indicate the corresponding error on the l a st dec imal. Based on 
photometric morphology, SDSS provides a separation between galaxy-like (3) an d star-like objects (6), see iLupton et al.l ( 1200 ih . *"' Frequency 
V = 1.4 GHz The column lists R2. If not available, Rl or Bl is given instead (see lMonet et al.|[2003l and references therein). 



6. Discussion 

6.1. An AGN origin of 1 FGL J0030. 7+0724 

The y-ray signal of IFGL J0030. 7+0724 can be explained 
by a conventional AGN. With respect to the unified 
scheme for the spectral energy distribution (SED) of AGN 
(namely F lat Spectrum Radio Quasars (FSRQs) and blazars), 
see, e.g., iDonato et al.l (1200 ih . the hard spectral index of 
1FGLJ0030.7+0724 (F ^ 1.7) is compatible with a high- 
energy-peaked blazar (HBL). Within the updated positional 
uncertainty of IFGL J0030.7+0724 derived from the 24- 
month data (Fig. [8]), the most likely radio counterpart is 
NVSS J003 119+072456 (/i.4ghz = (11-6 ± 0.6)mJy), which 
positionally coincides with the newly discovered hard X-ray 
source SWIFT J0031 19.8+072454 (F ^ 1.6). Note that corre- 
sponding to the notation of Table |2l the Swift source is flagged 
with a G in Fig. [8] The energy flux observed between 0.2 and 
2keV is ~ 2 X 10"'^ ergcm"^ s"' (Table|7]i. Additionally, an op- 
tical counterpart of the radio and X-ray source is listed in the 
SDSS catalogue (r = 17.4™), see Table [8] In Fig. |9]we show 
an empirical model for the average SED of HBLs, which is 
based on the bolometric luminosity distribution of FSRQs and 
blazars iFossati et al. 1997, 1998; Donato et al. 2001). The SED 
is normahsed to the radio flux of NVSS J003 119+072456 (at 5 
GHz). For comparison, the spectral measurements of the opti- 



cal and X-ray counterparts are presented as well. Within the ob- 
servational errors and assuming temporal variability, the y-ray 
spectrum of IFGL J0030. 7+0724 is consistent with the model 
prediction. Furthermore, the spectral index of the X-ray source 
agrees with an HBL, while its flux is fainter than predicted for an 
(average) HBL. This might be also explainable by temporal vari- 
ability (the radio. X-ray, and y-ray observations were not taken 
simultaneously), and blazars are well known to be variable in all 
wavelength bands, where the a mplitude of variability increases 
with energy (lUMch et al.lll997h . 



The other fainter objects in the uncertainty region (the 
radio source NVSS J003030+072132 and the two X-ray 
sources E and F, see Fig. [8]) are less likely to be associ- 
ated with IFGL J0030.7+0724, but cannot be excluded. For 
NVSS J003030+072132, no X-ray association was detected 
with Swift-XKY at the level of 2 x lO^'^^ergcm"^ s"'. No con- 
clusive optical counterpart is catalogued (above ~ 26"", see Sect. 
15.1.11 ). With respect to the comparatively high y-ray signal (cf., 
Fig-HJ' this source therefore fails to provide a convincing coun- 
terpart for IFGL J0030.7+0724. Similarly, the lacking radio de- 
tection as well as energy fluxes (~ 4x lO"''* erg cm"- s"'), which 
are much fainter than the HBL prediction, disfavour a coinci- 
dence of the X-ray sources E and F with IFGL J0030.7+0724. 
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Fig. 9. Energy spectra of IFGL J0030.7+0724 (solid red line) 
and the favoured radio (filled black triangle), optical (violet 
crosses, dereddened ugriz magnitudes), and X-ray (solid blue 
line) counterparts, together with the SED of an average HBL 
(solid black line). The SED was adapted from iDonato et al.L 
assuming the average redshift of known HBLs z — 0.25 
dPonato et alJ 1200 ih . and is normalised to the radio flux of 
NVSS J003 119+072456. The frequency-dependent energy flux 
v/v is given in the observer's frame. Note that the statistical er- 
rors of the radio and optical data points are too low to be re- 
solved in the figure. Statistical uncertainties of the X- and y-ray 
spectra are indicated by the corresponding s haded areas, which 
we derived with Eq. 1 in lAbdo et al.l (l2009bh . The filled red cir- 
cle indicates the catalogued high-energy flux from ferm/-LAT. 
Observations with ROSAT provide an upper limit on the X-ray 
flux at the nominal position of IFGL J0030.7+0724 between 0. 1 
and 2.4 keV, wh ich is depicte d by the green square (95% c.l., as- 
suming r = 2.0, lBormll20Tol) . 



6.2. A DM subhalo origin of 1FGL J0030.7+0724 

Without a clear indication for variability, it remains plausible 
that the y-ray emission of IFGL J0030. 7+0724 originates from 
a DM subhalo. The analysis of the arrival times of the source 
photons (Sect. 15.1.21 ) is consistent with a temporally constant 
source of moderate spatial extent. The reconstructed high-energy 
flux within the statistical errors is 0p(lO- 100 GeV) > 5 x 
10 " cm - s while the upper limit on the extent is 6^ < 0.7°, 
corresponding to 068 S 0.3°. As shown in Sect. |4] in realis- 
tic WIMP scenarios the high effective self-annihilation cross 
section required to explain the source with DM annihilation 
in a FHM subhalo is hardly compatible with current observa- 
tional constraints (see Fig. [3] and Table [1]). However, given the 
more realistic SHM, flux and extent of IFGL J0030.7+0724 are 
consistent with a subhalo of mass between lO*" and 10** Mq. 
Assuming a DM subhalo of 10*" M©, the resulting distance 
would be 2.4^Q ^' kpc, given the concentration scatter of the SHM 

model. For a WIMP of 500 GeV annihilating to bb, the required 
minimum effective enhancement is 7 for a high-concentrated 
SHM subhalo with a corresponding distance of 1 .7 kpc, while 
it increases to 31 for an average-concentrated subhalo with a 
corresponding distance of 2.4 kpc. Note that h{0.7°) x 1.4. 
An even lower boost factor is required for a lighter WIMP 
of 150GeV which predominantly annihilates to t^t : 3(13) 
for a high-concentrated (average-concentrated) SHM subhalo. 



Table 9. Fluxes above the energy thresholds of MAGIC and 
H.E.S.S., predicted by a DM scenario of IFGL J0030.7+0724. 



Flux prediction for MAGIC/H.E.S.S. [%Crab] 

150 GeV 500 GeV 1 TeV 
bh 0.3/ 10-' 0.6/0.05 

W*W- 0.5/0.01 0.8/0.2 

T+r- 0.7/- 3.1/1.1 



Notes. The fluxes are listed in percent ages of the^Crab Nebula's flux, 
0cmb(> 50 GeV) a; 1.6 X lO"** cm '^ s'' l|Alb ert et al. 2008a) and </>crab(> 
300 GeV) a 1.5 x lO"'" cm"^ s"' jAharonian e t aL 2Q06c). respectively. 
Effective cross sections required by the individual DM scenarios are 
discussed in Sect. 16.21 raised by a factor of about 2.3 for ra^ = 1 TeV. 



Further decrease of the necessary boost may be provided by sub- 
substructure and cuspier profiles (Sect.|4]i. 

In addition to theoretical uncertainties on halo properties and 
their expected scatter (Sect. ITTl i. observational uncertainties af- 
fect the distance and boost factor estimates. The uncertainties 
on the flux directly change the boost, while uncertainties on the 
most crucial measurement, the angular extent 0^, affect both the 
required boost and the distance estimate (Sect.[3]l. The discussed 
object IFGL J0030. 7+0724 serves as an appropriate benchmark, 
because the corresponding uncertainties are representative for a 
typical DM subhalo source. The observational uncertainties are 
of similar magnitude as the theoretical ones. 

6.3. Remarks and prospects for lACTs 

The 28-month data of Fermi-LAT contains no additional photon 
detected around the nominal position. This lowers the probabil- 
ity of steadiness to ~ 25% and may indicate variability, which 
supports a BL Lac scenario. Vice versa, such a behaviour would 
also be anticipated by a selection bias; If the true flux is lower 
than the value found in the discovery data set, the discovery con- 
dition would only allow for the detection of sources where the 
photon number has been fluctuating upwards. Poisson fluctua- 
tions of this faint signal could have accounted for a detection 
with the LAT even if the true flux had remained below the detec- 
tion sensitivity. 

It is instructive to note that with regard to a definite identi- 
fication of a counterpart (or ruling out a candidate) from obser- 
vations in other wavelength regimes the limiting factor is the 
accuracy of the ferm/-LAT source position {0(5'), cf.. Table 
|4]i and PSF. With just six detected photons, probably includ- 
ing one background photon, the source is close to the confu- 
sion limit. This situation can only be resolved by future instru- 
ments with much larger effective areas, such as the proposed 
CTA, which will probe deep into the expected population of sub- 
haloes. The much larger number of photons would help to infer 
significantly improved source positions. Furthermore, for detect- 
ing a spectral cut-off and in case of heavy WlMPs (;7i^ > 1 TeV), 
observations in the VHE range with lACTs are favoured. For 
the particular DM scenarios proposed for IFGL J0030.7+0724, 
fluxes anticipated in the energy ranges accessible for MAGIC 
and H.E.S.S., (/>(> 50 GeV) and (/>(> 300 GeV), respectively, 
are listed in Table |9] (given by E) oc Ny{> E)/Ny{lQ- 
100 GeV), see Eq. [8J. Additionally, flux estimates for WlMPs 
with nij^ - 1 TeV were derived. Note that the required effective 
cross sections (see Sect. 16.2b increase by a factor of 2.3, because 
<o-v)eff oc Vy(lO-lOOGeV)"'. Also note that fluxes expected 
for VERITAS are comparable to those for H.E.S.S. 
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The low energy threshold of MAGIC leads to comparatively 
high integrated VHE fluxes for < 1 TeV. The flux prediction 
for MAGIC is of 0(1%) of the Crab Nebula's for the favoured 
T^T-,m^ = 150 GeV and W+W-,m^ = 1 TeV model. With 
MAGIC, 50 hours of observation are necessary to detect this 
source with more than 5cr. For comparison, predicted fluxes for 
H.E.S.S. are not higher than 0.2% Crab for these models , which 
require s a few hundred hours of observation ( Aharonian et al.l 
l2006d) . We remark that advanced analysi s methods improve the 
sensit ivity of H.E.S.S. by a factor of 2 (Ide Naurois & Rollandl 
|2009|) . In the near future, an additional telescope (H.E.S.S.- 
II) will lower the energy threshold of H.E.S.S. to about 25- 
50 GeV. For the corresponding flux level of l%Crab, the re- 
quired observation time for H.E.S.S. -II and MAGIC will be sim- 
ilar. Furthermore, the planned CTA obse rvatory will be able to 
detec t such a source in about 50 hours (iThe CTA Consortium! 

l2oToh . 

7. Summary and conclusions 

Hierarchical structure formation predicts Milky Way-sized 
galaxies to host numerous DM subhaloes with masses rang- 
ing from 10'° down to a cut-off scale of lO"-' - 10"" Mq. 
Given standard WIMP scenarios, e.g., motivated by supersym- 
metry, we have demonstrated that DM subhaloes are detectable 
with the currently operating y-ray telescope Fenni-LALT. Based 
upon state-of-the-art models, detectable subhaloes would obser- 
vationally appear as faint high-energy y-ray sources between 
10 and 100 GeV with a flux at the sensitivity level of Fermi- 
LAT (~ 10"'° cm"^ s"' between 10 and 100 GeV for one year). 
The observable y-ray emission exhibits a moderate spatial ex- 
tent below ~ 0.5°. Subhaloes favoured for detection are mas- 
sive (10^ - 10^ Mq) at distances of (9(kpc), while low-mass 
subhaloes are not detectable. Within the intrinsic halo-to-halo 
scatter, only a moderate enhancement of the self-annihilation 
cross section preferred by standard cosmology, (crv)o = 3 x 
10"^^ cm"* s"', between 3 and 12 is necessary (dependent on the 
WIMP model), which is consistent with current observational 
constraints. Increasing sensitivity for a data-taking period of five 
years will allow us to resolve subhaloes requiring a cross sec- 
tion enhanced by a factor between 1.3 and 5. Additional sub- 
substructure within a subhalo may lower the required enhance- 
ment. Within statistics, one massive subhalo could be detectable 
with Fermi-V^NI in one year and might appear in the first-year 
catalogue (IFGL), assuming a subhalo population predicted by 
numerical A^-body simulations. Regarding the IFGL, the high- 
energy flux (10- 100 GeV) of subhalo candidates should be 
fainter than ~ 4 x 10"'° cm"^ s"' (for the WIMP scenarios con- 
sidered here). 

Intensive searches for subhaloes in the IFGL reveal twelve 
candidates, which are unassociated, non-variable, high-latitude 
sources detected above 10 GeV. The physical origin of the 
most promising object selected by lacking association, faintness, 
and spectral index, IFGL J0030. 7-1-0724, was investigated by 
analysing the 24-month data set of Fermi-LAT. With dedicated 
Swift-XRT observations (10.1 ks), seven X-ray sources were dis- 
covered around IFGL J0030. 7+0724. Located within the posi- 
tional uncertainty of the y-ray source, a radio source positionally 
coincident with a newly discovered X-ray source hints at a con- 
ventional HBL origin of IFGL J0030.7-H0724. However, owing 
to a large positional uncertainty and the lacking detection of vari- 
ability, the possibility of a dark nature remains. The measured 
high-energy flux and spatial extent of the source is compatible 
with a DM subhalo between lO*" and 10^ M© in a distance of 



about 2kpc, driven by a 500 (150) GeV WIMP self-annihilating 
to bb (t^t"). In this case, the required enhancement of (crv)o is 
7 (3) within the intrinsic scatter of the subhalo model, given a 
subhalo of 10'' Mq. 

Establishing the - probably more likely - HBL scenario of 
IFGL J0030. 7+0724 requires a significant detection of y-ray 
variability and a confirmation of the radio as well as X-ray coun- 
terparts. Vice versa, a steady y-ray flux with a spectral shape pre- 
dicted by self-annihilating WIMPs would hint at a DM nature of 
the object. This validates the necessity of additional intense and 
long multi-wavelength observations. In particular, lACTs offer 
a unique capability to reduce the positional uncertainty of faint 
LAT sources and to detect a spectral cut-off in the VHE range. A 
detection of the subhalo candidate IFGL J0030. 7+0724 may be 
possible with telescope systems like H.E.S.S. -II, MAGIC, and 
CTA within about 50 hours of observation. 

Our results encourage the search for more subhalo candi- 
dates in current and upcoming (very-)high-energy data releases. 
However, even in optimistic scenarios the expected number of 
LAT-detectable subhaloes is small. Furthermore, a longer expo- 
sure time - while certainly helpful with regard to the single can- 
didate discussed in this work - will not neccessarily remedy the 
general problem of the y-ray photon count that limits the po- 
sitional accuracy and therefore the chance of identifying coun- 
terparts. Given < 1 TeV, acquiring a sufficiently large num- 
ber of detections which may solve the subhalo problem requires 
higher sensitivity in the high-energy range. 
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Appendix A: Concentration of Aquarius subhaloes 

The Aquarius simulation provides results on the profile param- 
eters o f resolved subhaloes, taking tidal interaction into ac- 
count (iSpringel et a l. 2008a). These results are used to derive the 
distance-averaged virial concentration of subhaloes to confront 
it with the toy-model predictions used here. 



13 



H.-S. Zechlin et al.: Dark matter subhaloes as gamma-ray sources and candidates in the first Fermi-hAT catalogue 



Following up on Eq. [T] th e tidal concentration q = ^t/'"s 
is introduced (cf.. lAndoll2009l) . where Rt denotes the tidal and 
therefore physical radius of a subhalo. For an NFW-type mass 
density profile, q = exp[W(-e'") + a]- l,a= I + MJ{Anp^rl), 
where Wix) denotes Lambert's W-function and Mt the tidal sub- 
halo mass. In numerical simulations, the directly "observable" 
quantities of (sub)haloes are related to the dynamics of the halo 
system, including the maximum velocity Vmax and the distance 
'"max where Vmax is reached. To recover the canonical parame- 
ters and related to the density profile, we use approximate 
relations 2[y™ax/(//ow)]- - 5.80 x 10^ [Mt/(10'* Mo)]-^ '** 
and Mt ^ 3.37 x 10^ \ V^^J(\Okms- '^ )f^'^U^, fitting th e 
results of the simulation " (ISpringel et al.l l2008at lAndol l2009h . 
Given analytical relations betw een (rmax, Vmax) and (?"s,ps) for 
the NFW profile (e.g., Eq. 8 in iKuhlen et"al]|2008l) . this yields 
^ 0.094 [M,/(10^Mo)]"-3^kpc and pt\Mr) ^ 9.6 x 
Therefore, the tidal concentration 
cf'' is determined via a'^'^ ^ I + 0.66 [Mt/(10<' Mo)]""4, which 
is valid for masses above the resolution limit of the simulation, 
Mt > 3.2 X IO'^Mq. 

The virial concentration of Aquarius subhaloes is given by 

[3Mvir/(4;rAcPcrit)]'^V'"f''(^vir), where the charac- 
teristic radius as function of the virial subhalo mass is obtained 
from an empirical relation mapping Mvir to Mt. Based on cf and 
assuming the FHM virial concentration-to-mass relation (Eq.O, 
the relative tidal mass is Mt/Mvir ~ /[cf''(Mt)]//[c™(Mvir)], 
since M = Anp^rlf{c). Hereby, we have assumed that tidal ef- 
fects on inner subhalo parts are negligible: ps(Mvii)rs(Mvii)^ ~ 



lOVcrit [Mt/(10''Mo)]-o i«. 



C^?(Mvir) 



Ps(Mt)r,(Mt)^^ 



Mi-Myij; relation is 



With /(cf^) = 



7Aq _ 



1, the distance-averaged 



Mt(Mvir) ^ 



712.6 kpc 

47rpc: 



-3 



1.04 



/(c^™) 

^ vir ^ 



1.04 



(A.l) 



For massive subhaloes (> 10^ Mg), /(c^,';*^)"' "'* is well fit by 
a power law, f[c™^{M^ir)]-^°'^ ~ 034 [Mvh-/(10'' Mo)]°02, 
yielding MJM^i, ™ 0.23[Mvh/(10^ Mo)]^^^ for Mt > 3.2 x 
IO'^Mq. This reveals rf''(Mvir) ^ 0.054 [Mvi,-/(1O''M0)]"-4O kpc 
and, therefore, the distance-averaged virial concentration of sub- 
haloes 

0.07 

) ^ 46.8 1 i (a.2) 



10^ 



for Mvir 6 [1.5 X 10^;- 1O'°]M0. 

Appendix B: Moderately extended Fermi sources 

For y-ray catalogues such as IFGL, instrument data have been 
analysed assuming sources to be point-like. Given that de- 
tectable subhaloes would appear as moderately extended accord- 
ing to the PSF of Fenni-LAT (see Sect. 14X11 crpsp ~ 0.15° 
for E = lOGeV), we investigated the effect of the IFGL point- 
source-analysis framework on extended sources. 

To study the high-energy flux 0p(lO-lOO GeV) reconstructed 
by the point-source analysis for a given intrinsic (subhalo) extent 
0a, a Monte-Carlo (MC) simulation dedicated to the particular 
source IFGL J0030. 7+0724 was used. Based on the 11-month 
data set (see Sect. 15.1.21 for details), the celestial coordinates of 
each of the five source photons between 10 and 100 GeV were 
re-simulated. The intensity profile was assumed to follow the 
line-of-sight integral over the (squared) NFW profile of a sub- 
halo for the given 6*5 (peaking at the nominal source position). 




0.4 0.5 0.6 
es [deg] 

Fig.B.l. Average (solid lines) and rms (shaded areas) of 
0p(lO- lOOGeV) (upper panel), TS (middle panel), and the 
scaling li{6s) (lower panel) as function of the intrinsic angu- 
lar extent 9^. For each 6^, a sample of 500 simulations of the 
IFGL J0030.7-H0724 photon distribution between 10 and 100 
GeV assuming a DM subhalo intensity profile was analysed with 
the IFGL point-source-analysis framework. 



Other observational photon parameters, such as energy, inclina- 
tion, detection time, conversion type, and event class (see Table 
|6]l, were kept fixed. Subsequently, detectional influences were 
accounted for by smoothing with the PSF. For each 9^, 500 it- 
erations were analysed with the framework described in Sect. 
15.1.21 {gtfindsrc and gtUke) according to flux and significance 
{S ~ VTS, where TS denotes the test statistic of the analysis). 
All other sources within the ROI were kept fixed. The study is 
restricted to the signal-dominated regime chosen to be ^ 1° 
given the low background A^ba- Since a: 0.46°, this corre- 
sponds to ~ 3 cTpsp. Justified by the low background, all photons 
were treated as signal events. 

The 9s dependence of the sample-averaged reconstructed 
flux 0p(lO - lOOGeV) and corresponding test statistic TS is 
shown in the two upper panels of Fig. IB. II For large 9^, the prob- 
ability of photons to be located far away from their central po- 
sition increases. Therefore, both 0p and TS decrease because of 
a minor contribution of outer photons to the point-source region 
(defined by the PSF). For 9^ ~ \° , the average significance drops 
below the detection criterion (TS > 25). Note that TS > 25 still 
holds for about 35% of the simulated samples. 

In terms of Eq. [8] appropriate investigation of candidates 
provided by point-source catalogues is therefore admitted by a 
scaling h{0s), which allows us to map the catalogued flux 0p to 
the true flux (p of the entire source. The angular dependence of 
h is shown in the lower panel of Fig. IB. II Given - h(0s) (f>p, 
the factor was derived by defining h{Q°) = 1 . Conservatively, the 
complete MC sample was used to derive li(9s), including real- 
isations with TS < 25Q As expected in the signal-dominated 
regime, the increase of h with increasing 0^ is comparatively 
slight, while it is fairly linear in the background-dominated 
regime. Note again that this result holds for sources similar to 
IFGL J0030.7+0724 at high galactic latitudes only, while in 
general h = h{l, b, 0^,). 



Given the selection bias of the IFGL catalogue, TS > 25, a more 
stringent deduction of h{0s) should include realisations with TS > 25 
only. This lowers the effective scaling factor h(ds). 
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Vice versa. Fig. IB. II states a reasonable (but conservative) 
value of the sensitivity of Fenni-LAT for hard sources of similar 
type: 0p(lO- lOOGeV) * lO^'^cm^^g-i ^^jg t hat this value 
is sim ilar to the point-source sensitivity stated in lAtwood et al.l 
(l2009l) . 



Appendix C: Subhalo-induced diffuse flux 

In the following, the diffuse flux of th e subh alo population 
is derived using a prescription by lAndol (l2009h . which is ex- 
tended to include the probability distribution of the concentra- 
tion parameter c (see Eq. |4]i. Numerical A^-body simulations 
have demonstrated that the diff'erential subhalo number density 
d«sh = M{r, M) dM follows a power-law in subhalo mass M. 
Following standard assumptions, the number density N(r,M) 
factorises such that N{r, M) oc n^^iir) ■ M where a - 1.9 and 
r is the distance to the host's centre. In simulations, the spa- 
tial density distribution n^'nir) is consistently found to be "anti- 
biased" and, e.g., «sh('") <^ PEin('") ( Springel et al 2 008a) . where 
PEinir) denotes the Einasto profile (lEinastdll965b 



P¥.m{r) oc exp 



— I -1 

r-2 



(C.l) 



For a Milky Way-sized halo, the best-fit parameters for the sub- 
haloes' spatial distribution pv\„{r ) have been found to be - 
0.68 and r_ 2 ^ 0.81 c^^rf^ dSoringel et al.l l2008al) . where 
c^* ^ 15 ( ICatena&Ullioll20T0l) . Using N{r,M) normalised 
to represent a probability density function in M, the differential 
density is 



= «sh('")- 



dM 



M 



(C.2) 



where M^-^m ^ M^^^ are the minimum and maximum mass of 
Galactic subhaloes, respectively. The normalisation of the sub- 
halo number density rishir) is chosen such that the fraction of the 
host's mass distributed in subhaloes /sh = M^^^jM^^ =15% for 
the cut-o ff scale M„in = 10''' M a, where M^* = (1 .49 + 0. 17) x 
10' - Mq (ICatena & Unioll2010l) . The chos en value of f^h is con- 
sistent with recent estimates /sh - 10-50% ('Diemand et a l.l2005t 
Piema nd & Moore 2009; Diemand et al. 2008b; Spring el et all 
2008bh . The total mass contained in subhaloes is given by 



■ An 



I Arr^ \ 
Jo Jm,, 



dMM 



d« ,h(r, M) 
dM 



(C.3) 



Solving for an Einasto-type profile and a yields 



3/ffE-l 

— I r 



3 2 (R] 



2-a 



(a - 1)(A- 



2-a . 



1) 



exp 



as o'E V '■-2 
O'E I r-2 



-1 



(C.4) 



where r(a, x) is the lower incomplete gamma function and A - 

Mmax/Mmin. 

The minimum mass M^nin of subhaloes is governed by the 
details of kinetic decoupling of WlMPs in the ea r ly Universe 
(jBerezinskv et al.1 12003L 120061; iGreen et al.1 120051; iBringmannI 
I2009h . Depending upon the mass and composition of, e.g., the 
neutralino, a wide range of minimal subhalo masses has been 
considered in the literature, namely Mmin e [10'"; lO'-'lM©. 
Here, two benchmark cases are considered for Mmin, i.e., 10''° 



and 10 ^ Mq, bracket i ng the 500 GeV neutralino scenario dis- 
cussed bvlBring mannI (l2009l) . The upper mass limit was fixed to 
Mn,ax = lO'^M^* ^ 10"^ Mq. Results do not depend on the 
exact value chosen for Mmax- 

Using X.{M,D) (Sect. IZTI and I2.2l i. the average specific in- 
tensity from a subhalo population with extended, isotropic emis- 
sivity profiles is given towards a galactic direction h by 



«)) =JdM J ds- 



ds 

Mmm s,(X(M,j)) 



dnf,ii{r(s,h),M) 



dM 



dcP{c,c) 



ZAM,c) 

4-K 



, (C.5) 



assuming that the spatial extent of each subhalo is much smaller 
than the scale on which the subhalo distribution changes sig- 
nificantly. The total photon rate is given by Eq. |6] with the 
substitution Ny — > EdNy/dE. Furthermore, £.v{M) is required 
to be one-to-one. The galactocentric radius corresponding to 
the position sh is r{s,^) - (Rq + - 27?oicosi/')'^^, where 
i/r denotes the angle between h and (cos(t/') = {h,Ro}). 
Subhaloes bright enough to be detected as individual sources are 
not considered to contribute to the diffuse emission. Therefore, 
the lower limit of the line-of-sight integral is set by the detection 
criterion X. > 47rij0sens, where (p^^ns denotes the flux sensitivity 
for a detection in one year with ferm/-LAT, see Sect. |4] Since 
Ro «: Ryii, the upper bound of the ^-integral imax(«) ~ R^'^ ■ The 
SHM photon rate is a function of both M and s, and therefore s, 
also depends slightly on s. Conservatively, s,{M, s) = s,(M, s), 
s = n^'^ revealing a lower bound on i,. 

vir ' ^ 
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